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The Hubbard model for electrons with orbital degeneracy is shown to have an underlying 5'(7d(4) 
symmetry of spin-orbital double. A hidden charge S(7c(4) symmetry is exposed and an extended 
Lieb-Mattis transformation which maps these two symmetries into each other is given. On the basis 
of elementary degenerate perturbative theory, it is shown that the system with strong repulsive 
coupUng is equivalent to a SO(6) Heisenberg magnet at half-filling and a SU(4) one at quarter-filling. 
The band is half-filled at all temperature for ^ = 3(7 /2. The features of ground state and low-lying 
excitations in one dimension are indicated according to exact solutions. Various possibilities of 
symmetry breaking are also given. 

PACS number(s): 71.10.Fd, 71.10.-w, 75.10.Jm, 71.30. -Fh 



I. INTRODUCTION 

There has been much interests in the studies on corre- 
lated electrons in the presence of orbital degree of free- 
dom [|l| because the orbital degree of freedom plays an 
important role in understanding the phenomena, such as 
metal-insulator transitions, high-temperature supercon- 
ductivity and colossal magneto-resistance. The orbital 
degree of freedom is relevant to many transitional metal 
oxides It may be also relevant to some Cfjo ma- 

terials [7| and samples of artificial quantum dot arrays 
For a theoretical understanding of the observed un- 
usual properties, a SU(4) theory describing spin systems 
with orbital degeneracy was proposed [p|. There were 
also numerical |l^ and perturbative jll| studies of 1- 
dimensional models for these systems. The ground-state 
phase diagrams for the system with a symmetry break- 
ing of S';7(4) ^ S';7(2) X SU{2) were discussed ]ll|,@. 
Experimentally, the phase separation |]l3| was observed. 
Due to the rapid developments in experiments where the 
metal ions has orbital degeneracy in addition to spin de- 
generacy, a theoretical study of such system by taking 
account of the kinetic terms caused by nearest neighbor 
hopping becomes indispensable. 

In this paper we study a Hubbard-type model for elec- 
trons with orbital degeneracy. In Sec. g, we show that 
the model has an underlying SU(4) symmetry of spin- 
orbital double. The spin and orbital operators are re- 
lated to the SU(4) generators, which will be helpful for 
further studies on the magnetization. A hidden charge 
SU(4) structure is exposed in Sec. |l|. An extended 
Lieb-Mattis transformation which maps those two SU(4) 
symmetries into each other is also presented. From a 
basic relation derived from particle-hole transformation, 
we show that the band is half-filled at all temperature 
when the chemical potential equals 3C//2. In Sec. IV, 
two kinds of "partially negative ?7" models are intro- 
duced and analyzed according to the strategy of p4| . 
Three basic excitation modes are show to exist in the 
spin-orbital sector. In Sec. ^ with the help of the par- 



tially attractive models, we study the repulsive large U 
model on the basis of elementary degenerate perturbative 
theory. It is shown that the effective Hamiltonian with 
strong repulsive coupling at half-filling is equivalent to 
the Hamiltonian of SO (6) Heisenberg model, and that at 
quarter-filling is equivalent to the one of SU(4) Heisen- 
berg model. In Sec. ^ several poss ibili ties of SU(4) 
symmetry breaking are given. In Sec. VII, we summer- 



ize the main results and discuss the agreement with exact 
solutions in one dimension by describing the features of 
ground state and low- lying excitations. 



II. UNDERLYING SU(4) SYMMETRY 

We consider electrons with doubly orbital degeneracy. 
The spin components are denoted by up (f) and down 
(I), the orbital components by top and bottom. The 
four possible states of electrons are 
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We use 1, 2, 3, and 4 to label these states from now on. 
Let us consider the Hamiltonian of electrons with two- 
fold orbital-degeneracy on a lattice 

H = -t Ct{x)Ca{x') + Y.^--''^-^''>-'^''y (2) 

<x,x'> 

where x's identify the lattice site, a, a' = 1,2,3,4 spec- 
ify the spin and orbital as defined in the above. The 
C+(x) creates a fermion of state \a) located at x site and 
na{x) is the corresponding number operator. Eq. (^) is 
the Hamiltonian for four-component systems, and there 
were various discussions on multi-component Hubbard 
model in one dimension |jl^,|l6|. We remark that the 
four-component Hamiltonian can also describe either a 
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spin-3/2 system jlj or a toy model of proton and neutron 
system with on-site strong interaction. In the terminol- 
ogy of group theory, the former forms a high-dimensional 
(here it is 4 dimensional) representation of Ai Lie alge- 
bra, while the later forms the fundamental representation 
of D2 Lie algebra. Whence, the physics that Eq. (||) de- 
scribes will be precise only when the representation space 
for the internal degree of freedom is specified. It refers 
to the spin and orbital in our present discussion. 

We can verify that the Hamiltonian (^) with Uaa' = U 
commutes with following 15 operators, 

= \ - C++i(x)C„,+i(a;)], 

X 

Ea^ = X! '^m(2;)C,„+l(x), 

X 

E-a.^ = {EaJ\ m= 1,2,3, (3) 
and additionally, _Eq,j+q2 = [Ea-i^,Ea2], Ea^+as — 

(-^01-1-02)^? Q2^Ct3 — (-^Ct2 +Ct3 ) ^ 7 and £/_Q^_Q2— as — 

+0,3+0.3)^, here am's stand for the simple roots of 
A3 Lie algebra. These operators precisely obey the com- 
mutation relations of At, Lie algebra ||l8|. Whence the 
Hamiltonian (|^) with Uaa' = U is invariant under any 
global SU(4) rotation. We denote this underlying sym- 
metry of spin-orbital double by SUdi'i)- There are several 
equivalent ways to write out the generators of the Lie al- 
gebra. We adopt the Chevalley basis because physical 
quantities can be conveniently evaluated in this basis. 
Because the A3 Lie algebra is of rank three, i.e., three 
generators in its Cartan sub-algebra, there are three con- 
served quantum numbers that label the eigenstates. In 
terms of these 0„i's that commute to each other, the z- 
components of both total spin S^^^ and total orbital T^^^ 
are give by 

T/„, = 01+202+03. (4) 

They are useful for evaluating magnetizations once the 
ground state in the presence of magnetic field is solved. 



III. HIDDEN CHARGE SU(4) SYMMETRY 

In addition to the above symmetry, one may easily 
think of the U{1) charge invariance jl^. Moreover, there 
exists a larger hidden symmetry in present model on a 
bipartite lattice, we call charge S'C/c(4) symmetry. Their 
Chevalley bases are given by 

X 



X 

F^,=Y,e'^-^C+{x)C+{x), m=l,2,3, (5) 

X 

where x = (tt, tt, ...), Qm's are generators of Cartan sub- 
algebra of the A3 Lie algebra. The other generators are 
given by standard relations that we demonstrated previ- 
ously when observing the underlying SUd{^) symmetry. 
The charge SUc{'l) symmetry is not only valid for the 
Hamiltonian (||) at high energy scale t ^ U but also 
valid for a kind of on-site coupling which we are going to 
show. Considering U13 = U24 — —U while Uab = U for 
the other subscripts and taking account of the chemical 
potential term in the Hamiltonian, we obtain 

[H',F^J = (-l)'"(2/i-C/)i^„,„, 

H' = H - fiJ^M^)- (6) 

x.a 

The commutators between H' and Qm's always vanish. 
Apparently, the model has a charge 5?7c(4) symmetry 
when fi = U/2. The mentioned requirement for the sign 
of the on-site coupling constants is unnecessary for the 
traditional Hubbard model which has a hidden charge 
SU{2) symmetry because there is only one constant 
for coupling of spin up and spin down. 

Eq. (^) implies that the raising operators F^^ of the 
charge SUc{^) create some pair of electrons to a given 
state. Precisely, F^^ or ^0,3 creates a double occupancy of 
spin singlet, however, F-a^ creates a double occupancy of 
spin triplet. These operators provide mappings between 
states in distinct sectors of different electron numbers. 

There exists an extended Lieb-Mattis transformation: 

C,{x)^e'^-^Ct{x), z = 2,4, 
C,{x)^C,{x), J = 1,3, (7) 

which maps S'C/c(4) into 5L/d(4) and vice versa. The ap- 
plication of particle-hole transformation gives rise to a 
basic relation |l^ for bipartite lattice: 

E{Na,U)^E{L-Na,U) + i{N ^2L)U, (8) 

where L is the total number of lattice sites and N the to- 
tal number of electrons. Using this relation we can derive 
a relation for the thermal average: 

< N >^,T= 4L- < N >3U-t.,T, (9) 

where N — J2x a''^<^i^)- ^ result, the band is half- 
filled at all temperature when /i = 3C//2. 

IV. PARTIALLY ATTRACTIVE MODELS 

The magnitude and sign of the on-site coupling may 
vary from system to system. For isotropic pure attractive 
coupling Uab = U < 0, the unperturbed ground state has 
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iV/4 of the sites occupied by "quaternarys" . The ground 
state is degenerate when 7V/4 < L because the energy 
3NU/2 is independent of which sites are occupied. As a 
generahzation of Cooper pair, the quaternary might have 
abundant physical meanings. It can form a SU(4) singlet 
for Uab = U . It can also form either two pairs being spin 
singlet but orbital triplet or that being spin triplet but 
orbital singlet, or form a "resonance" state being alterna- 
tions of them depending on the symmetries remained in 
a general Uab- For example, a triplet-pairing supper con- 
ductivity model was discussed by a particular choice 

of Uab- 

We have particular interests in two kinds of partially 
attractive on-site couphngs: (i) U13 = C/24 = —U while 
Uab = U for the others; (ii) U12 = U23 = U13 — —U while 
Ui4 = U24 = U34 = U. Here U < for both cases. Let 
us consider these two cases respectively. 

In case (i), the local favorite states in energy are a 
quaternary and four types of pairs. The quaternary is an 
instantaneous state which may separate into two pairs 
randomly because there is no difference in energy be- 
tween them. The hoping terms in the Hamiltonian split 
the degeneracy to form a band of charge-density wave 
states where the quaternary and pairs move from site 
to site. As is known in the absence of orbital degree of 
freedom fl^ that the spin-density excitation turns over 
a spin to break a pair at a cost in energy. The spin- 
orbital-density excitation, however, turns over either a 
spin or a orbital which results in three basic processes. 
In addition to the process of breaking a pair at the cost 
\U\ in energy, one process transmits a favorite pair into 
an unfavorite pair at the cost 2\U\ in energy. The process 
of breaking a quaternary into a trinity and a single costs 
\U\. Breaking a quaternary into two unfavorite pairs will 
cost 4|[/|, which is not a basic process because it can be 
represented by two processes of the second type. All the 
other complicated processes can always be represented 
as a composition of those three basic processes. Thus we 
believe that there are three elementary quasi-particles in- 
volved in the excitations in the spin-orbital sector. 

For the case (ii), the favorite states in energy are two 
types of trinitys and three types of pairs. The charge- 
density wave states arising from the hopping terms are 
of the movement of the trinitys and pairs from site to 
site. The spin-orbital-density excitation that turns over 
either a spin or a orbital involves in three basic processes 
which cost \U\, 2\U\ and zero in energy. Therefore, there 
are three elementary quasi-particles, in which a gap-less 
node is expected to exist. 

With the help of the above discussions we are now in 
the position to employ elementary degenerate perturba- 
tion theory for more quantitative formulations. It is suf- 
ficient for calculating the low temperature properties to 
consider the lowest band only. For those two types of par- 
tially attractive on-site couplings that were discussed pre- 
viously, all electrons remain in either the favorite pairs, 
trinitys, or quaternary in the lowest band. Because it 
breaks pairs, trinitys, or quaternary, the perturbation 



part (the hopping term) has vanishing first order matrix 
elements. Thus it must be calculated to second order 
by considering the virtual transitions into the next band. 
After some algebra one gets, 



ie 



^o)ag = J2<9\Hf\g'>ag 

a' 

_ < g\Ht\^lj > 

dq — 

£ - £0 



(10) 



where \g > denote various degenerate states of the unper- 
turbed ground states, i.e., Hij\g >= £o|(7 >. Ht stands 
for the hopping term and Hu for the interaction term 
of Eq. (^. After the similar calculation as in ||l^, we 
obtain the effective Hamiltonians that will be helpful for 
studying repulsive model. 



V. EFFECTIVE MODELS OF STRONG 
REPULSIVE COUPLING 

We now study the repulsive model Uab = C/ > 0. 
First we consider the half- filled band [N = 2L). For 
the ground state in this case, every sites are doubly occu- 
pied by electrons. The excitation of charge-density waves 
brings about inevitably a triple occupancy of site at a 
cost of \U\ in energy at least. By making a particle-hole 
canonical transformation: 



2,4, 
1,3, 



Cj{x) ^ Cj{x), j 
the Hamiltonian becomes 

H=t {-irc:ix)Caix') 

<x,x'> 

+ UY^ {-ir+'''na{x)na' (x) + Va{na{x) + i), (11) 



where = [1 — 3(— l)"]C//2. Clearly, the repulsive on- 
site coupling becomes the partially attractive case (i) 
which we already discussed. In the unperturbed states, 



the sites originally occupied by 
whereas those occupied by 



i 



i 



T 



T 



> become empty, 
> are replaced by 



T i 



T i 



>. The other four kinds of double occu- 



pancies exchanged their positions. The degenerate per- 
turbation theory can now be used in the same way as 
it was done in the partially attractive models that we 
considered before. After reversing the canonical trans- 
formation, we obtain the effective Hamiltonian as follows 



h(x, x') 



\U\ 

E 

ran 



E 



h{x, x') 



g 



— a. (j^ )■ (12) 
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where A denotes the set of roots of the A3 Lie alge- 
bra, Ea{x) and Ora{x) are generators of the Lie al- 
gebra given by Ea^{x) = C+{x)Cm+i{x), Om{x) = 
C+{x)C,n{x) - C:^^-^{x)C„i+i{x), m = 1, 2, 3. These op- 
erators generate the whole states from the highest weight 

T i 



state 



> in the following way, 



T i 
1^- 



> 



i T 



> 



> 



E-a, 



\ 



> 



E-ai 



i 



T i 



> 



where the x in Ea{x) is omitted. Those states form the 
vector representation of A3 (isomorphic to D3) Lie alge- 
bra. Whence it is a SO (6) Heisenberg model. 

As to the case of quarter-filling (N=L) with repulsive 
on-site coupling, the ground state has one electron on 
each site. The excitation of charge-density waves brings 
about double occupancy of sites indispensably and costs 
energy of \U\ at least. It is again interesting to make a 
canonical transformation 

Ciix) ^ C+ix), Qix) ^ Ciix), I = 1,2,3. 

Consequently, the sites at first occupied by | | > are 
now empty whereas the other sites are doubly occupied, 
i.e., by either \ ]_ J >, | A >, or | ^ J >. The re- 
pulsive on-site coupling becomes the partially attractive 
case (ii) which we discussed previously. Again, the degen- 
erate perturbation theory can be applied as before. The 
effective Hamiltonian is obtained as a SU(4) Heisenberg 
model. The local states form the spinor representation of 
A3 Lie algebra, and their relations in the weight diagram 
were given [pi by: 



E^r 



I > I i > E-^' I Y > I I 



E^r 



> 



The evaluation of correlation function in the strong cou- 
pling limit will be considerably simple by means of the 
corresponding effective Hamiltonians. 



VI. SYMMETRY BREAKING 

The underlying SU(i(4) symmetry is fulfilled for the 
isotropic on-side coupling Uaa' — U. There is no phase 
separation between spin wave and orbital wave at the 
SU(4) point. A complete separation is expected to occur 
after the anisotropic on-side couplings in the spin-orbital 
configuration are introduced. Since the 'diagonal' cou- 
pling Uaa has no physical contribution due to fermionic 
wave function vanishes when two electrons at the same 
site being the same SU(4) state (i.e., a=a'), we can intro- 
duce six parameters Vat (a < b) to break down the SU(4) 
symmetry, i.e., Uaa' = U + Vaa'- 

It is not difficulty to find the possible symmetry break- 
ings by calculating the commutation relations between 
the Hamiltonian and the SU(4) generators. There are 
two ways to break the SU(4) down to the SU(3)xU(l). 
For V12 ^ V23, vi4 = V24 = V34 (or V12 = V13 = V14, 

V23 — V2A ~ 1^34), the residue symmetry SU(3)xU(l) 
is generated by {Oi, O2, £^±ai , ): £^±021 O3} 

(or by {Oi,02,03,£:±c2,£^±(a2+Q3)£^±a3})- This is 
a two-parameter hierarchy. For a three-parameter 
hierarchy, ■i;i3 = W23 = ui4 = U24, the 
residue symmetry is SU(2)xSU(2)xU(l) generated by 
{Oi,E±a^,OzjE±a3i02}. When V12 = W13 and V2a = 
1^34, the residue symmetry becomes SU(2) xU(l) xU(l) 
with {O2J -E'±Q2 7 Oi, O3} as its generators. This is obvi- 
ously a four-parameter hierarchy. Furthermore, if either 
V12 7^ t'13 or 1)24 7^ W34 the previous SU(2) xU(l) xU(l) 
will be broken to the U(l) xU(l) xU(l) generated by 
{01,02,03}. 

As the SU(4) Lie algebra is of rank 3 (there are three 
generators in its Cartan subalgebra) the Zeemann-like 
interactions with external fields reads 

Hz = y^femO„i(a:), 



where m — 1,2,3. In general, it breaks the 
symmetry down to the minimum residue symmetry 
U(l)xU(l)xU(l). However, if h2 = = 0, it regains 
a SU(2)xU(l)xU(l) generated by {d, 02, O3, £;±a3}. 
Similarly another SU(2) xU(l) xU(l) generated by 
{Oi,E±at02, O3} regains for hi /i2 = 

Since it is the SU(4) singlet, the ground state is invari- 
ant under any SU(4) rotations. Except for the singlet- 
excitation and the pure charge-excitation states which 
arc still of invariant under SU(4) rotation, the other 
multiplet-excitation states vary. For the system with 
N — ^n, there is always an axis in flavor space along 
which the multiplet-excitation states are invariant under 
a rotation of 27r. 



VII. DISCUSSIONS 

We have shown both the underlying and hidden sym- 
metries of Hubbard model with orbital degeneracy. We 
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derived the effective Hamiltonian in strong repulsive cou- 
pling for both half-filled and quarter-filled band. The 
band is half-filled at all temperature if the chemical po- 
tential is 3J7/2. It is shown that there are three elemen- 
tary modes involved in the excitations in the spin-orbital 
sector. 

We did not specify the dimensions in the above dis- 
cussions on the symmetries of the system and its low- 
temperature effective Hamiltonians. It is interesting to 
consider one dimensional case |Q because the exact so- 
lution in one dimension always provides non-perturbative 
features. The exact solution can be obtained by means 
of Bethe-Yang ansatz similar to ref. if states with 
site occupation of more than two are excluded ^4 \ . Al- 
though the multi-component generalizations of Hubbard 
model in one dimension were explored by several authors 
in various aspects, the relationship between the 
enlarged internal degree of freedom and concrete physics 
was not clearly exhibited. It is convenient to take ther- 
modynamics limit to the Bethe ansatz equation by in- 
troducing density distributions of the quasi-momentum 
k and that of three rapidities for the spin-orbital double 
(we call SU(4) flavor). The ground state being of real 
roots and no holes is a SU(4) singlet, accordingly, both 
spin and orbital are in "anti-ferromagnetic" states. The 
features of low-lying excitations are studied by consider- 
ing the contributions of holes and complex 2-strings. An 
important consequence is that the excitations in charge 
sector and flavor sector are separated. In charge sector 
the elementary modes are holon and anti-holon (parti- 
cle), and the real excitations are gap-less particle-holon 
and gapful holon-holon excitations etc. There are three 
types (in agreement with the above general analyses in 
any dimension) of flavorons as the elementary excitation 
modes in flavor sector, namely, two quadruplets carrying 
spin 1/2 and orbital 1/2 that form respectively the fun- 
damental representation and the conjugate representa- 
tion of SU(4); a hexaplet carrying either spin 1 or orbital 
1 which forms a 6-dimensional representation moreover. 
These flavorons compound to constitute real excitations. 
The details of one dimensional model are given in another 
paper [^ . 

We also analyzed the possible symmetry breakings 
caused either by extending the model to anisotropic on- 
side coupling or by introducing external field. Our anal- 
yses based on the Hamiltonian structure will be helpful 
for further discussions on the phase diagram by means of 
numerical DMRG. 
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